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Abstract: By analyzing the large-angle Bhabha scattering events e*e ^ — (y)e*e^ and diphoton 
events e*e^ — yy for the data sets collected at center-of-mass (c.m.) energies between 2.2324 and 
4.5900 GeV (131 energy points in total) with the upgraded Beijing Spectrometer (BESIII) at the 
Beijing Electron-Positron Collider (BEPCII), the integrated luminosities have been measured at the 
different c.m. energies, individually. The results are the important inputs for R, value and J/w 


resonance parameter measurements. 


Key words: luminosity, Bhabha, diphoton, R value 


PACS numbers: 13.66.De, 13.66.Jn 
I. INTRODUCTION 


Hadron production in e*e^ annihilation is one of the 
most valuable testing grounds for Quantum Chromody- 
namics (QCD), and is an important input for precision 
tests of the Standard Model (SM). The R value, which is 
defined as the lowest level hadronic cross section normal- 
ized by the theoretical +7 production cross section in 
e*e- annihilation, is an indispensable input for the de- 
termination of the non-perturbative hadronic contribu- 
tion to the electromagnetic coupling constant evaluated 
at the Z pole (a(M2)) [1, 2], and the anomalous mag- 
netic moment a, = (g — 2)/2 of the muon [3]. The domi- 
nant uncertainties in both a( M2) and a, measurements 
are due to the effects of hadronic vacuum polarization, 
which cannot be reliably calculated in the low energy 
region. Instead, with the application of dispersion rela- 
tions, experimentally measured R values can determine 
the effect of vacuum polarization. 


In experiment, the R value is determined by 


obs bk 
R = Nees — Nyad (1) 
otp L- ena EE. (144) 


where N$5$ is the number of observed hadronic events, 


NS is the number of background events, £ is the inte- 
grated luminosity, &n4q is the detection efficiency for the 
hadron event selection, £15 is the trigger efficiency, 1+6 
is the initial state radiation (ISR) correction factor, and 
c^, is the Born cross section of ete~ — ^j^. There- 
fore, the measurement of integrated luminosity plays an 


important role in the R value measurement. 


Quantum electrodynamics (QED) processes can usu- 
ally be used to determine the integrated luminosity due 
to larger production rates, simpler final state topologies 
and more accurate cross section calculation in theory rel- 
ative to the other processes. The integrated luminosity 
is measured by 


Nil — NS 


po Nein = Nar _ (2) 
rig 
OQED * EQED * EQED 
where Ng is the number of the QED events observed 


in the experimental data, NEED is the number of back- 


ground events, cQgp is the cross section of the selected 


QED process, corn is the detection efficiency and QD 


is the trigger efficiency. 


In this paper, we present the measurements of lu- 
monisities of the R scan data samples taken at BESIII 
from 2012 to 2014. The measurements are performed 
by analyzing two QED processes e*e^ — (y)ete~ and 
e*e- — yy. For energy points near the J/«» resonance, 
only the e*e^ — yy process is used, because J/y) — 
(y)e*e- events can not be distinguished from e*e^ > 
(y)e*e- events experimentally. 


II. DETECTOR 


BEPCII [4] is a double-ring e*e^ collider designed to 
provide a peak luminosity of 10?? cm~?s~! at the center- 
of-mass (c.m.) energy (ys) of 3770 MeV. The BESIII [4] 
detector has a geometrical acceptance of 93% of 4x and 
has four main detector sub-components: (1) A small-cell, 
helium-based (60% He, 40% C3Hg) main drift chamber 
(MDC) with 43 layers providing an average single-hit 
resolution of 135 um, and charged-particle momentum 
resolution in a 1 T magnetic field of 0.596 at 1 GeV/c. 
(2) An electromagnetic calorimeter (EMC) consisting of 
6240 CsI(TI) crystals in cylindrical structure arranged in 
a barrel and two end-caps. The energy resolution at 1.0 
GeV/c is 2.5% (5%) in the barrel (endcaps), and the po- 
sition resolution is 6 mm (9 mm) in the barrel (endcaps). 
(3) A time-of-flight (TOF) system for particle identifica- 
tion composed of a barrel part made of two layers with 
88 pieces of 5 cm thick, 2.4 m long plastic scintillators 
in each layer, and two endcaps with 96 fan-shaped, 5 cm 
thick, plastic scintillators in each endcap. The time res- 
olution of 80 ps (110 ps) for barrel (endcap) prodvides 
2e K/m separation for momenta up to ~ 1.0 GeV/c. (4) 
A muon system (MUC) consisted of 1000 m? of resistive 
plate chambers in nine (eight) layers of barrel (endcap) 
provides 2 cm position resolution. 


IIl. DATA SAMPLE AND MONTE CARLO 
SIMULATION 


The measurements of luminosities are performed for 
131 data samples, including 4 energy points at 2.2324, 
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2.4000, 2.8000, 3.4000 GeV taken at the 2012 run, 104 
energy points from 3.8500 to 4.5900 GeV taken at the 
2013-2014 runs, 15 energy points near the J/~ produc- 
tion threshold, 4 energy points during the 7 mass mea- 
surement and 4 energy points for charmonium studies. 


The ete~ — (y)e*e^, yy and (y)ut u^ events are sim- 
ulated with the generator Babayaga v3.5 [5]. The back- 
ground process of e*e^ — r*T- is generated with the 
KKMC [6], while the e*e^ — hadrons and ete~ — ete- 
+ X (X can be hadrons or leptons) events are generated 
with LUARLW [7] and BesTwogam [8], respectively. 


IV. ANALYSIS 


The e*e^ — (y)ete~ events are required to have two 
good charged tracks with opposite charge. Each charged 
track is required to be within +10 cm of the interaction 
point in the beam direction and 1 cm in the plane per- 
pendicular to the beam. In addition, the charged tracks 
are required to be within |cos6| < 0.8, where 0 is the 
polar angle, in the MDC. Without applying further par- 
ticle identification, the tracks are assigned as electron 
and positron depending on their charges. The deposited 
energies of electron and positron (+) in the EMC are 
required to be larger than 0.65 x Ebeai to suppress back- 
grounds, where Eyeam is the beam energy. To make sure 
the the selected charged tracks are back-to-back in the 
c.m. system, |A@+| = |01 + 05 — 180°| < 10.0? and 
[A^ó.-| = ||ó1 — ¢2| — 180?| < 5.0? are required, where 
61/2 and $,;5 are the polar and azimuthal angles of the 
two charged tracks, respectively. Figure 1 shows the com- 
parisons of the momentum and polar angle distributions 
of electron and positron between experimental data and 
Monte Carlo (MC) simulation at \/s = 2.2324 GeV, the 
good agreements are observed. 


To select e*e^ — yy events, the number of good 
charged tracks is required to be zero. Two neutral clus- 
ters are required to have a polar angle | cos0| < 0.8 with 
the deposited energy E, satisfied 0.7 < E,/Epeam < 
1.16. The two selected photon candidates are further 
required to be back to back by applying the requirement 
|Ag,| = |ó51—6423| < 2.5°, where $,1/2 are the azimuthal 
anlge of the photons. Figure 2 shows the comparisons of 
the enegy deposition, polar angle and Ao, distributions 
of two selected photons between experimental data and 
MC simulation at ys = 2.2324 GeV. 


The numbers of observed QED events, Neb: are ob- 
tained by event-counting after applying the event selec- 
tion requirements on experimental data at different c.m. 
energies, individually. The detection efficiencies of sig- 
nals, €qgEp, are obtained by analyzing the corresponding 
signal MC events as done in data analysis. The cross sec- 
tions of selected QED processes are calculated with the 
Babayaga v3.5 generator and the trigger efficiencies are 
quoted from Ref. [9]. 


To estimate the numbers of background events, Nor: 


two different methods are applied for e*e^ — (y)e*e^ 
and e*e^ — yy processes, individually. In e*e^ > 
(y)e*e^ process, the numbers of background events are 
estimated by performing the same requirements on the 
background MC samples, which yields a background level 
of 1075 after normalization. In e+e~ — yy process, the 
background level is relatively large due to the hadronic 
process contamination. The normalized numbers of back- 
ground events from e*e^ — yy are estimated from the 
Aó., sideband region, defined as 2.5? < |Aó,| < 5.0. 
The distributions of the A@, sideband is supposed to be 
flat by analyzing the background MC samples. 


Table I shows input numbers used to calculate the lu- 
minosities at \/s = 2.2324 and 3.0969 GeV. 


TABLE I. Summaries of the input numbers in luminosity cal- 
culation at ys = 2.2324 and 3.0969 GeV. 


JVs | QED | Nib [NOH] conv aa £ 
(GeV) | process (nb) | (96) | (96) (pb) 
(y)e*e |728522| 8 1476.5|18.74| 100 | 2.645 
2.2324 yy 86974 | 1138 | 70.26 |46.50| 100 | 2.627 
3.0969 yy 36083 | 1062 | 36.59 |46.25| 100 | 2.069 


V. SYSTEMATIC UNCERTAINTY 


The main systematic uncertainties of the integrated 
luminosity are originated from the uncertainties related 
to the requirements on the kinematic variables, tracking 
efficiency, cluster reconstruction efficiency, c.m. energy, 
MC statistics, background estimation, trigger efficiency 
and generators. 


For the systematic uncertainty from requirements on 
each kinematic variable, we re-measure the luminosity 
by altering the required values, i.e., | cos0| < 0.8, | A0, 
< 10°, [Age] < 5°, |Ad,| < 2.5°, E;« / Eseam > 0.65 
and 0.7 < Ey/Ebeam < 1.16, individually. The resultant 
differences of measured luminosity with respective to the 
nominal value are taken as the systematic uncertainty. 


To study the uncertatinty of tracking efficiency, a 
Bhabha event sample is selected with only EMC informa- 
tion [10]. The candidate events are selected by requiring 
the two clusters registered in the EMC with the deposited 
energy larger than 0.65 x Epeam and lied within the polar 
angle | cos0| < 0.8, corresponding to the angular cover- 
age of the barrel EMC. Since the two clusters originated 
from e* in the e*e^ — (y)ete~ candidate events are 
bent in the magnetic field, the two shower clusters in the 
xy-plane of the EMC are not back-to-back. Aó.« is re- 
quired to be in the range of [—40°, —5?] or [5?, 40?] to 
remove the e*e^ — yy events. We further apply the 
MDC information on the selected candidates, and the 
ratio of survived events is regarded as the tracking effi- 
ciency. The average difference on the tracing efficiency 
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FIG. 1. The distributions of momentum (up plots) and polar angle cos 0 (down plots) for electron (left) and positron (right) at 
V/s = 2.2324 GeV. Dots with error bars are experimental data and red histograms are signal MC simulation. The MC entries 


are normalized to the experimental data. 


between data and signal MC simulation, 0.4196, is taken 
as the systematic uncertainty. 


The systematic uncertainty due to the cluster recon- 
struction efficiency in the EMC is determined to be 0.0596 
for e^ by comparing the cluster reconstruction efficien- 
cies between data and signal MC (both for et and e7). 
Since high-energy y and e^ behave in good approxima- 
tion in the EMC, the value of 0.0596 is also taken as the 
systematic uncertainty due to the cluster reconstruction 
efficiency in the EMC for a single y. 


The uncertainty of c.m. energy is estimated to be 2 
MeV [11]. For each energy point, an alternative MC 
simulation sample of 1 million events with a c.m. en- 
ergy of 2 MeV above the nominal value are generated 
to re-estimate the detection efficiency, the results differ- 
ence is regarded as the systematic uncertainty from c.m. 
energy. 


The uncertainty of MC statistics is 0.1796 for the e* e^ 
— (w)e*e- process and 0.1596 for the e*e^ — yy pro- 
cess, which is estimated by 


(3) 


where N is the number of signal MC events, and e is the 
detection efficiency. 


The rate of background events in the selected e*e^ — 
(y)e*e^ candidate events is very small (1075). There- 
fore, the uncertainty due to background contamination is 
neglected. For e* e^ — yy events, the rate of background 
events is the normalized number of selected background 
events in the sideband region divided by the number of 
signal events, which are (1.53+0.03)% and (1.31+0.04)% 
for experimental data and the MC simulation, respec- 
tively. Therefore, the difference 0.2396 is taken as uncer- 
tainty from background contamination. 


The trigger efficiencies for barrel e*e^ — («)e*e^ 
events and e*e^ — yy events are 100% with an uncer- 
tainty of less than 0.196 [9]. 


'The uncertainty due to the Babayaga generator v3.5 is 
0.596 for ete~ — (y)e*e-, while 1.096 for e* e^ > yy [5]. 


Systematic uncertainties at ys = 2.2324 GeV for e*e^ 
— (vy)e*e- and e*e^ > yy are listed in Table II. As- 
suming all sources of systematic uncertainties are uncor- 
related, the total uncertainty is calculated to be 0.7% for 
e*e- — (y)e*e- and 1.1% for e*e^ — yy by adding all 
the contributions in quadrature. The uncertainties re- 
lated with the tracking efficiency, cluster reconstruction 
efficiency, trigger efficiency and generators are common 
between the different c.m. energy points, while others 
are c.m. energy dependent and are determined for the 
different c.m. energy points, individually. 
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Deposited energy distributions of the most energetic y (upper left), the secondary most energetic y (upper right), 


cos (bottom left) and Ag (bottem right) at ys = 2.2324 GeV. Dots with error bars are experimental data and red histograms 
are signal MC simulation. The MC entries are normalized to the experimental data. The discrepancy in the deposited energy 
distributions is due to the imperfect simulation of energy correction deposited in TOF. However, it will not affect the efficiency 
since loose requirements on these variables are applied. The uneven distribution of cosÓ is due to the structure of crystals in 


the EMC. 


TABLE II. Summary of systematic uncertainties at ys = 


2.2324 GeV. 


Source 
|cos0| < 0.8 
[Aba] « 10° 
[Ages] < 5° 
[Aġ,| < 2.5? 

E. / Ebeam > 0.65 
E.- / Eseam > 0.65 
0.7 < E^, | Beam < 1.16 
Tracking efficiency 
Cluster reconstruction 
Beam energy 
MC statistics 
Background estimation 
Trigger efficiency 
Generator 
Total 


VI. 


By using the QED p 


et 


TABLE III: The summaries of measured integrated luminosities from the 
two QED processes. The first uncertainty is statistical and the second 


SUMMARY 


sured for 131 data samples with c.m. energy between 
2.2324 and 4.5900 GeV. The precision of integrated lu- 
minosity is around 0.7% for ete~ — (y)ete~, while 
around 1.1% for ete~ — yy. The total luminosity is 
1036.3 pb-!, and the luminosities at the individual c.m. 
energy point are summarized in Table III. The ratio of 
the measured luminosity from two process is illustrated 
in Fig. 3. The ratios are closed to 1 within the uncer- 
tainties, which indicates the results from the two mea- 
surements are consistent well with each other. For each 
energy point out of the J/w resonance region, the lu- 
minosity measured by ete~ — (y)ete~ is more precise 
and thus is recommended. For energy points around J/w 
(from 3.0930 to 3.1200 GeV), only the luminosities mea- 
sured by e* e^ — yy are obtained. The measured results 
are the important inputs for the physics studies, e.g., R 
value measurement and J/1 resonance parameter mea- 
surement. 


one is systematic. 


rocesses e*e^ — (»y)e*e- and 


€ — yy, the integrated luminosities have been mea- 


ete — (y)e*e 


(pb 


7 


TABLE III: (continued) The summaries of measured integrated TAREE III: (continued) The summaries of measured integrated lumi- 
nosities from the two QED processes. nosities from the two QED processes. 
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Vs (GeVy[e > (Mere (pb) ere 2 17 (pb-7] Vs (GeV) [ete — (ere (pb 
2.8000 3.753+0.008+0.026 3.766+0.014+0.042 4.0500 
3.0500 14.893+0.030+0.103 14.919+0.029+0.158 4.0550 
3.0600 15.040+0.030+0.131 15.060-0.029-+0.158 4.0600 
3.0800 31.019+0.060+0.189 30.942+0.044+0.338 4.0650 
3.0830 4.740+0.011+0.029 4.769+0.017+0.052 4.0700 
3.0900 15.709+0.031+0.099 15.558+0.030+0.162 4.0800 
3.0930 14.910-+0.030+0.157 4.0900 
3.0943 2.143+0.011+0.023 4.1000 
3.0952 1.816+0.010+0.019 4.1100 
3.0958 2.135+0.011+0.023 4.1200 
3.0969 2.069+0.011+0.024 4.1300 
3.0982 2.203+0.011+0.023 4.1400 
3.0990 0.756+0.007+0.008 4.1450 
3.1015 1.612+0.010+0.018 4.1500 
3.1055 2.106+0.011+0.022 4.1600 
3.1120 1.720+0.010+0.019 4.1700 
3.1200 1.264+0.009+0.013 4.1800 
3.4000 1.754+0.012+0.020 4.1900 
3.5000 3.643+0.017+0.040 4.1950 
3.5424 8.711+0.027+0.098 4.2000 
3.5538 5.593+0.021+0.059 4.2030 
3.5611 3.894+0.018+0.043 4.2060 
3.6002 9.620+0.028+0.108 4.2100 
3.6500 48.618+0.065+0.538 4.2150 
3.6710 4.603+0.020+0.052 4.2200 
3.8500 7.962+0.028+0.088 4.2250 
3.8900 7.799+0.028+0.087 4.2300 
3.8950 7.626+0.027+0.085 4.2350 
3.9000 7.631+0.027+0.085 4.2400 
3.9050 7.625+0.027+0.085 4.2430 
3.9100 7.267+0.027+0.082 4.2450 
3.9150 '1.5332:0.027-E0.088 4.2480 
3.9200 6.903+0.026+0.076 4.2500 
3.9250 6.76340.026+0.075 4.2550 
3.9300 6.825+0.026+0.076 4.2600 
3.9350 7.144+0.027+0.079 4.2650 
3.9400 7.256+0.027+0.082 4.2700 
3.9450 7.608+0.028-+0.086 4.2750 
3.9500 7.739+0.028+0.086 4.2800 
3.9550 8.141+0.029+0.090 4.2850 
3.9600 8.548+0.029+0.095 4.2900 
3.9650 7.770+0.028+0.086 4.3000 
3.9700 7.368+0.028+0.082 4.3100 
3.9750 8.050+0.029+0.089 4.3200 
3.9800 7.808+0.028+0.087 4.3300 
3.9850 7.992+0.029+0.089 4.3400 
3.9900 8.104+0.029+0.091 4.3500 
3.9950 7.984+0.028+0.084 4.3600 
4.0000 7.805+0.028+0.088 4.3700 
4.0050 7.567+0.028+0.085 4.3800 
4.0100 7.164+0.027+0.079 4.3900 
4.0120 6.951+0.027+0.079 4.3950 
4.0140 6.716+0.027+0.075 4.4000 
4.0160 6.582+0.026+0.074 4.4100 
4.0180 6.996+0.027+0.078 4.4200 
4.0200 6.735+0.027+0.075 4.4250 
4.0250 6.58340.026+0.073 4.4300 
4.0300 16.526+0.042+0.187 4.4400 
4.0350 6.687+0.027+0.074 4.4500 
4.0400 6.640+0.027+0.073 4.4600 
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TABLE III: (continued) The summaries of measured integrated lumi- 


nosities from the two QED processes. 
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FIG. 3. The ratios of luminosities measured by e*e^ — 
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(y)e*e^ and ete” — yy. The major plot is for the data 
samples with c.m. energy larger than 3.8500 GeV, while oth- 
ers are shown in the insert plot, the two methods give fully 
compatible results within the quoted uncertainties. 
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